Theory of Resonant Inelastic X-ray Scattering by Collective Magnetic Excitations 
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I present a tractable theory for the Resonant Inelastic X-ray Scattering (RIXS) of magnons. The 
low-energy transition operator is written as a product of local spin operators and fundamental x-ray 
absorption spectral functions. This leads to simple selection rules. The scattering cross section 
linear (quadratic) in spin operators is proportional to the fundamental magnetic circular (linear) 
dichroic spectral function. RIXS is a novel tool to measure magnetic quasi particles (magnons) and 
the incoherent spectral weight, as well as multiple magnons up to very high energy losses, in small 
samples, thin films and multilayers, complementary to Neutron scattering. 
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Magnetic excitations in solids are traditionally inves- 
tigated with inelastic neutron scattering. This technique 
has led to a better understanding of magnetism in gen- 
eral and shaped the thinking of physicists in terms of 
magnons as the magnetic quasi-particles present within 
a solid. It is important that the interaction of neutrons 
with the magnetic moment is well understood and can 
be approximated by a function linear in spin operators. 

One might expect that x-rays are not suited for the 
study of magnetic excitations, as photons, contrary to 
neutrons, do not carry a magnetic spin and therefore can- 
not directly excite single magnetic excitations. Recent 
resonant inelastic x-ray scattering (RIXS) experiments 
see clear dispersing magnetic excitations • More high 
quality data is expected due to the experimental improve- 
ment in both brilliance and resolution [5, 6]. For neutron 
scattering the magnetic dipole interaction between the 
neutron and electron spin is well understood. For RIXS 
the effective magnetic interaction is highly non-trivial. 
The optical dipole transition does not introduce spin- 
flip excitations. The magnetic excitations are a result of 
core-hole spin-orbit coupling and many-body interactions 
active in the resonant intermediate state [7Hl4j. 

Based on local cluster calculations De Groot et al. [7] 
realized that magnetic excitations are in principle possi- 
ble. Luo et al. showed that, within the fast collision 
approximation, it is possible to describe the RIXS spectra 
in terms of a low-energy scattering operator, which can 
excite single magnons. At the Cu £2,3 edge this approx- 
imation allows one to describe the magnetic excitations 
including their momentum and polarization dependence 
[jj [Toj| . In general however the RIXS spectra show not 
only a strong momentum and polarization dependence, 
but also a strong resonant energy dependence [l|, LUL( . In 
order to include this, one has to take into account the in- 
teractions in the intermediate state, without neglecting 
the momentum dependence of the final state. 

Within this letter I present an effective low-energy op- 
erator for RIXS, focused on magnetic excitations. The 
idea of the derivation, is not to make approximations to 
the Hamiltonian, or to the Green's functions used to eval- 



uate the RIXS spectral function. Using symmetry argu- 
ments an exact solution for magnetic excitations can be 
obtained by factorizing the RIXS transition operator into 
symmetric polynomials of spin operators and fundamen- 
tal x-ray absorption spectral functions. The fundamental 
spectral functions describe the dependence of the x-ray 
absorption on the magnetic moment m. These are the 
isotropic, magnetic circular dichroic (XMCD) and mag- 
netic linear dichroic (XMLD) spectral functions [l5| . 

I will start the derivation with the Kramers-Heisenberg 
formula for the double differential cross section: 
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where i (f) labels the set of initial (final) states, Ei is 
the energy of state H is the Hamiltonian, and T is 
the life-time broadening, which, if all states are included 
in the Hamiltonian goes to zero. In practice one only 
keeps the relevant states in the basis set of H. The in- 
fluence of the neglected states is included by a possible 
state-dependent non-zero value of T. The photon energy 
and the polarization of the incoming (outgoing) photon 
is given by u^( ) and s^ \. T e is the optical transition 
operator given by the perturbation of the photon field 
on the system: T e = p • A, with A the photon vec- 
tor field and p the electron momentum operator. The 
state \i) represents the ground-state, for example a mag- 
netically ordered Neel state. The operator T creates a 
core-hole and adds a conduction electron to the system. 
The state given by T\i) is in general not an eigen-state 
of H and it will fluctuate in time. The Green's func- 
tion (uji + Ei + iT/2 — H)~ x describes this propagation. 
A possible propagation path that would add a magnetic 
excitation to the system could be the flip of the core spin 
due to the core-hole spin-orbit coupling. Many other ex- 
citations are possible. After the Green's function has 
acted on the state T\i) the operator fills the core hole 
and at the same time annihilates a valence electron. This 
entire process might leave the system in an excited state. 
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Using S(lj + Ei - E f ) oc 3[lim r ^ + (u + E { - E f 
zr/2) -1 ] one can introduce the following notation: 
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The operator i?^ £o is the low-energy effective scatter- 
ing operator. The operator both creates and annihi- 
lates a core hole, leaving the system in a low-energy 
excited state. The uji dependence of R £ J^° describes 
the true resonant nature of RIXS. The Green's function 
(u + Ei + iT/2 - H)- 1 in Eq. [2] describes how these 
low-energy excitations, for example a spin- flip, propagate 
through the sample. 

It is useful to factorize the momentum dependence of 
R s j £ ° . The photon field, A, can be expanded around 
each atom at site rj. One should take into account that 
A changes phase when one changes lattice position. For 
photon fields described by a plane wave the phase shift 
is given by a simple exponential: 



p-A = J2 elk ' rj Pr A 3 
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where the sum is over all atoms, k is the photon wave 
vector, rj is the position of atom j and pj (Aj) is the mo- 
mentum operator (photon vector field) expanded around 
atom j. 

Tj y£ (rt ) creates (annihilates) a core hole at site j 
(j f ). For core level spectroscopy one can assume that the 
core hole does not hop from one site to another. The site 
where the core hole is created is thus the same as the site 
where the core hole is annihilated. Inserting this in the 
definition of R £ J^° yields: 
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with q = ki — k Q . 

Before factorizing the effective scattering operator 
R^f q into spin operators and fundamental x-ray absorp- 
tion spectral functions I will address a few general prop- 
erties of Rfj £ q - It is interesting to note that the con- 
ductivity tensor at x-ray energies is defined such that: 
e* • <j • a = (i\Rfj^q\i)' The absorption for a given po- 
larization is proportional to — ^[e* • a • e] at q = 0. In 



order to calculate the inelastic scattering between state 
\i) and state |/) one needs to evaluate (f\Rfj £ q\i)- Using 
that2(/|^|i) = + + 

i) - (1 + i)((i\RZ} i e £\i) + {f\R £ J £o q \f)) one realizes that the 
RIXS transition probability for making excitations from 
state \i) to state |/) is determined by x-ray absorption 
spectral functions of the initial state, the final state and 
linear combinations of these two states. For magnetic 
excitations \i) could be a state with spin down, |/) a 
state with spin up. The real and complex linear com- 
binations are states with rotated spin directions. The 
fundamental spectral functions describe the dependence 
of the x-ray absorption on the magnetization direction 
of a given system. For magnetic excitations the RIXS 
transition probability between state \i) and state |/) is 
thus given by some linear combination of fundamental 
spectral functions. Below I will derive this dependence 
exactly with the use of symmetry arguments. 

For magnetic excitations one can restrict the basis set 
of R £ J £ q to the ground-state and all states with all pos- 
sible spin excitations from the ground-state. Within this 
basis one can use the theorem of operator equivalence as 
introduced by Stevens [16] for the description of crystal- 
fields in Rare-Earth compounds. Any operator acting in 
spin-space only can be written as a polynomial of spin 
operators. R £ J £ ° = P(S), where P remains to be deter- 
mined. 

The terms in P are restricted by symmetry relations. 
I will first assume spherical symmetry and later branch 
down to the possible crystal symmetries. In general R £ J £o j 
must be a scalar. In spherical symmetry R £ J £ ° can only 
depend on the angle between and e*. Operators of the 
form el • e<, e£ x e % • 5, (eg • 5)(e< • S) + (e< • S)(e* • S) and 
(eg • €i)(S • S) are allowed. 

The pre-factors of these different symmetry-allowed 
terms can be obtained using the notion that the RIXS 
spectral function for zero energy loss should reproduce 
the elastic scattering. The equation for the elastic scat- 
tering follows from the Kramers-Heisenberg equation by 
restricting the possible excitations to the ground-state: 
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At the same time, for elastic scattering one can take 
advantage of the optical theorem which states that the 
elastic scattering is related to the x-ray absorption spec- 
tral function |17j |. The elastically scattered signal is pro- 



portional to 



J , confirming the relation between 



and the conductivity tensor. Hannon et al [18[ 
related the elastic intensity to the magnetization direc- 
tion; this can be used to determine the pre-factors for 
the inelastic scattering. The elastic scattering is given 
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by 0,0: 



(i\RX-\i} = a^6 i -6* + a^e i xe* -m j 
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a^dei ■ mj )(e* ■ rrij) 



where the numerical value of the complex isotropic 

spectral function, cr^ and cr^ are the complex, fun- 
damental, XMCD and XMLD spectral functions respec- 
tively and rrtj is a unit vector in the direction of magneti- 
zation. The operator R £ J £ ° is given in terms of symmetric 
spin operators as: 
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where Sj is the spin operator acting at site j and s the 
expectation value (Sj) = s(s + 1). For real crystals the 
local point-group has to be included. This leads to a 
branching of the fundamental spectral functions [15| and 
polynomials of higher orders in S, up to an order of 2s 
hm. Eq. [8] is the main result of the present letter and an 
exact form of the RIXS transition operator for magnetic 
excitations truncated to single site transitions. It leads 
to beautifully simple selection rules: One can measure 
single spin-flip transitions with cross polarized light for 
spins in the plane of the polarizations. The operator 
quadratic in S can lead to single (AS = 1) and double 
(AS = 2) spin-flip transitions, as well as contributions 
to the elastic line. For 5=1/2 systems is zero by 
symmetry. 

In order to clarify the method I present RIXS spectra 
at the ^2,3 edge of samples containing Cu 2+ in sym- 
metry and Ni 2+ in Oh symmetry. I will use linearized 
spin-wave theory for a Heisenberg-model on a ID chain, 
a 2D square and a 3D cubic lattice with nearest neigh- 
bor interactions. Linear spin- wave theory assumes local 



ordered moments. For perfect one dimensional samples 
the ground-state is not ordered and for S = 1/2 systems 
the magnetic excitations are two spinon excitations [2o| . 
It is often more transparent to write R^j as the inner- 
product of two polarization vectors and a 3 x 3 tensor. 
In such a matrix notation, R^j for Cu 2+ , a S = 1/2 
system, in tetragonal symmetry becomes: 
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The fundamental spectra of Cu 2+ determining the 
RIXS transition at the £2,3 edge can be seen in Fig. la. 
They show a striking feature; the a^ } and a^J spectra 

a ig 

are zero. This is a direct consequence of having only one 
hole with x 2 -y 2 symmetry in which one can not excite a 
p electron with z polarized light. With RIXS one can not 
measure spin-flip excitations for fully aligned spins with 
S z = ±1/2 oriented perpendicularly to the d x 2_ y 2 or- 
bital [9]. As shown recently, direct spin- flip scattering is 
allowed in all other situations [10| . For real systems mag- 
netic excitations will always be measurable as S = 1/2 
fully aligned spins never exist in a solid. 

Knowing the transition probability one can calculate 
the RIXS spectral function, using for example linearized 
spin- wave theory [20]. In Fig. lc-e I show the spec- 
tral function as one would measure with RIXS in the 
cuprates. One finds that the intensity goes to zero at 
T and diverges at the antiferromagnetic Bragg peak. 
This is a well known behaviour for collective excitations 
[2q( | . The calculations for a two dimensional square lat- 
tice are in good agreement with resent measurements on 
La 2 Cu0 4 i- 

For Ni 2+ in cubic symmetry (a S = 1 system) the 
resonant inelastic scattering transition matrix becomes: 
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The numerical values of a at the £2,3 edge for Ni 2+ 
are plotted in Fig. lb. There are four fundamental spec- 
tra as the XMLD (<J^) spectra branch to two different 

(2) (2) 

representations (dej and &t 2 )• ^ n or der to separate the 
different contributions one can use different scattering 
geometries or different resonant energies. One should re- 
alize that also interference terms between the different 



channels have to be considered. The fundamental x-ray 
absorption spectra are well known and can be calculated 
with quite good accuracy, simplifying the disentangle- 
ment of the different scattering channels. 

In Fig. lc-h, I show the one (panels c-e) and two (pan- 
els f-h) magnon intensity as one could measure with RIXS 
in Ni 2+ samples. For real geometries one would mea- 
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FIG. 1: (color online) Left: Fundamental x-ray absorption spectra that enter into the RIXS transition operator as energy 
dependent complex matrix elements calculated for (a) Cu 2+ and (b) Ni 2+ . Right: The Cu 2+ and Ni 2+ one magnon (c,d,e) and 
Ni 2+ two magnon (f,g,h) RIXS spectral function, calculated using linear spin- wave theory for a ID chain (c,f), a 2D square 
(d,g) and a 3D cubic (e,h) Heisenberg model in energy loss units of z S J (number of neighbors x spin x exchange constant). 



sure a linear combination of the two spectral functions 
depending on resonance energy and polarization as dis- 
cussed above. For the single magnon spectral function 
only a single peak is seen at every momentum. For the 
two magnon spectral function a band is observed with 
the requirement that the sum of the momenta for the 
two magnons is equal to the transferred q. One therefore 
might expect not to see any q dependence. The single 
magnon intensity diverges at the momentum correspond- 
ing to the antiferomagnetic order. This leads to an ob- 
servable dispersion for the two magnon spectral function 
in one dimensional samples (Fig. If). For a three dimen- 
sional sample, the main peak of the two magnon spec- 
tra shows hardly any detectable momentum dependence 
(Fig. lh), in good agreement with recent measurements 
on NiO [2l|. The one magnon branch does disperse. 

In conclusion, I have presented a tractable theory for 
the calculation of RIXS spectral functions. It fully incor- 
porates the intermediate state Hamiltonian including the 
core-hole spin orbit coupling and the multiplet features 
as well as the momentum dependence of the low-energy 
excited states. This is achieved by factorizing the RIXS 
spectral function into two parts for which different ap- 
proximations are used. The low-energy Green's function 
can be calculated within any of the standard approxi- 
mations available to describe spin-waves or other exci- 
tations at hand. The transition operator is written as 
local operators multiplied by fundamental x-ray absorp- 
tion spectra. The transition operator has a simple form 
that allows one to determine selection rules depending 
on the polarization of the incoming and outgoing light. 
RIXS measurements can provide us with detailed infor- 
mation on both the magnons and the incoherent spectral 
weight, which contains information on the momentum 
dependent interactions of magnons with themselves and 
other degrees of freedom in the system. The RIXS cross 
section at the TM £2,3 edge is so large that the tech- 
nique can be used to measure small samples, thin-films 
or multi-layers making it possible to measure magnetic 



excitations in a whole new class of materials. 
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